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D. Valenti*,^ L. Magazzu^,^’^ P. Caldara,^ and B. Spagnolo 

^ Dipartimento di Fisica e Chimica, Group of Interdisciplinary Physics, 

Universitd di Palermo and CNISM, Unitd di Palermo, 

Viale delle Scienze, Edificio 18, 1-90128 Palermo, Italy 
^Radiophysics Department, Lobachevsky State University of Nizhni Novgorod, 

23 Gagarin Avenue, Nizhni Novgorod 603950, Russia 
(Dated: March 11, 2015) 

Normally, quantum fluctuations enhance the escape from metastable states in the presence of 
dissipation. Here we show that dissipation can enhance the stability of a quantum metastable 
system, consisting of a particle moving in a strongly asymmetric double well potential, interacting 
with a thermal bath. We find that the escape time from the metastable state has a nonmonotonic 
behavior versus the system-bath coupling and the temperature, producing a stabilizing effect. 
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Introduction. — Recently, the role of dissipation on 
the dynamics of quantum systems has been the subject 
of renewed interest [T]. 

The presence of a dissipative environment indeed influ¬ 
ences significantly the escape from a quantum metastable 
state. This is a general problem, of interest in many ar¬ 
eas of physics, whenever a sudden change in the state of 
a system occurs on timescales small with respect to the 
typical times of the systems dynamics. 

The archetypical model describing the escape process 
is that of a particle subject to a cubic or asymmetric 
bistable potential and linearly coupled to a heat bath of 
quantum harmonic oscillators mi. In such a system 
the decay from the metastable state occurs on timescales 
that depend on the friction and temperature. Various 
physical systems such as magnetization in solid state sys¬ 
tems H, proton transfer in chemical reactions d and 
superconducting devices [7] can be described within this 
framework. 

Calculations of the decay rates, using a cubic poten¬ 
tial, have been performed in Refs. di using functional 
integral techniques. In Ref. Hj starting with the par¬ 
ticle at the bottom of the metastable well, it has been 
shown that the decay rate decreases monotonically as 
the damping increases and grows with the bath temper¬ 
ature. Similarly, by using a master equation technique, 
a monotonic increase of the escape rate, with respect to 
the temperature, is found in Ref. m for a Gaussian wave 
packet initially in the metastable well of a biased quartic 
potential. 

Stabilization of a quantum metastable state by an ex¬ 
ternal time-periodic driving, in absence of environment, 
was obtained in Ref. m- Moreover, suppression of 
activated escape from a metastable state by increasing 
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the temperature was found in a time-periodically driven 
quantum dissipative system [12]. 

Common wisdom is that environmental fluctuations al¬ 
ways enhance the escape from a quantum metastable 
state. A critical issue and of great importance is: 
can the dissipation enhance the stability of a quantum 
metastable state? 

To answer this question we follow the time evolution of 
the populations of spatially localized states in a strongly 
asymmetric bistable system, starting from a nonequilib¬ 
rium initial condition. This choice allows us to observe 
how, increasing the damping, the relaxation process to¬ 
wards the stable well goes from a population transfer in 
which the metastable well is temporarily populated, to 
a mechanism of direct transfer to the stable state. This 
stabilization effect is related to that due to the suppres¬ 
sion of tunneling by dissipation in quantum regime mm- 
As a result we find that dissipation can enhance the sta¬ 
bility of the quantum metastable state. Indeed, we ob¬ 
serve that the escape dynamics is characterized by a non¬ 
monotonic behavior, with a maximum, as a function of 
the damping strength: there is an optimal value of the 
damping strength which maximizes the escape time, pro¬ 
ducing a stabilizing effect in the quantum system. This 
result, which resembles the phenomenon known, in the 
classical context, as noise enhanced stability (NES) of 
metastable states uni [H], sheds new light on the role 
of the environmental fluctuations in stabilizing quantum 
metastable systems. We also find that the behavior of the 
escape time versus the temperature is nonmonotonic with 
a minimum. Therefore as the temperature increases, an 
enhancement of the escape time is observed, increasing 
the stability of the quantum metastable state. 

Model. — We consider a quantum particle of effective 
mass M in a double well potential (see Fig. [^. The 
system’s bare Hamiltonian is Hq = p^/2M -\-V{q), where 
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Here wq is the natural oscillation frequency around the 
minima, AC/ is the barrier height and e the asymmetry 
parameter. In this work e is large enough to mimic the 
cubic potential, which is the archetypal model potential 
for metastable systems. 



FIG. 1: (Color online) Potential profile V{q) (Eq. Q) for 
AU — lAhuio and e = MTwJq. Horizontal lines; the first 

6 energy levels. Vertical lines: the position eigenvalues in the 
DVR. The dashed curve is the initial probability distribution 
|v['(a;, 0)p. For the tunneling splitting we have Ai? 4,3 = E 4 — 
E 3 = 0.2ftu)o while E 2 — Ei — 0.985fia;o- 

The environment is a thermal bath of N independent 
harmonic oscillators with position coordinates Xj and 
coupled with the particle through the linear interaction 
term 

In our study it is assumed that, in the continuum 
limit N —>■ 00 , the bath has the Ohmic spectral den¬ 
sity J(a;) = M^oj with a cut-off at a frequency much 
larger than any other involved in the model. The damp¬ 
ing constant 7 is a measure of the overall particle-bath 
coupling strength and has the same role as the classical 
damping in the quantum Langevin equation associated 
to the present problem [i] . 

In the quantum regime, given the particle’s initial 
preparation, the bath temperature is such that the dy¬ 
namics is practically confined among the first 6 levels of 
the potential shown in Fig. In this reduced Hilbert 
space, performing a suitable transformation, we pass 
to the discrete variable representation (DVR) [TS], in 
which the particle’s reduced dynamics is described in 
terms of the localized basis of the position eigenstates 
{ki), • ■ ■, ke)}, where q\q^) = qi\qi). 

Analytical method. — We assume a factorized initial 
condition, with the bath in the thermal state psifi) = 
I Zb- The particle’s reduced density operator in 
the DVR is given by 
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where Pfj.y{t) = {q^j.\p{t)\qu) {p.i' = l,---, 6 ) and 

K(qfj^,q^,t;qa,qp,0) is given by the double dissipative 


path integral 

rVqit) r V*q'{t)A[q]A*[q']EFv[q.q']. (3) 
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Here A[q\ is the amplitude associated with the path q{t) 
of the bare system. In the DVR a path is a step-like 
function with transitions among the positions qi. The 
effect exerted by the bath on the quantum mechani¬ 
cal amplitude associated to a path {q{t),q/{t)) is con¬ 
densed in the Feynman-Vernon (FV) influence functional 

J^Fv[q,q'] HU- 

This approach is non-perturbative in the system-bath 
coupling, and is thus suited for dealing with the strong 
coupling regime. Nevertheless, the FV influence func¬ 
tional makes the path integral intractable as it introduces 
time nonlocal interactions between the paths q{t) and 
qf{t), through the bath correlation function Q{t). 

The nonlocal part of the interactions cancels out in the 
limit in which the bath correlation function Q(t) is linear 
in t, i.e. in the long time limit t 3> h/ksT (see section IV 
and appendix B in Ref ini)- If the temperature is suffi¬ 
ciently high, Q{t) can be taken in the linearized form at 
all times, which amounts to perform the so called general¬ 
ized non-interacting cluster approximation (gNICA |17j . 
the multi-level version of the NIBA for the spin-boson 
model mm- By comparing the transition probabilities 
per unit time among the \qi)'’s with ksT/h, we obtain the 
limit T > O.I/kuo/fcB as a rough estimate for the validity 
of the gNICA for our system. 

Within the gNICA the double path integral of Eq. (§ 
assumes a factorized form in the Laplace space, allow¬ 
ing for the derivation of a generalized master equation 
(GME). If / 9 ( 0 ) is diagonal in the position representation, 
the GME reads m 
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Ptj,^i{t)='^ dt''H^^{t-t')pi,^{t'). (4) 

^= 1-^0 

The elements of H are taken to the second order in the 
transition amplitudes per unit time A^- = {qi\HQ\qj)/H 
and at all orders in the system-bath coupling. They dis¬ 
play a cutoff of the type exp(—yt x const.). At strong 
damping (as in our case, see below), in the short time in¬ 
terval in which is substantially different from zero, 
Piiu{t) is practically constant. This allows us to write 
the following rate equation as the Markov approximated 
version of Eq. @ 
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where = /“ dT'Hfj.^ir). 

The smallest, in absolute value, of the nonzero eigen¬ 
values of the rate matrix F determines the largest time- 
scale of the dynamics, that is the quantum relaxation 
time Trelax m- 
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Escape time. — In the following we focus on the par¬ 
ticle’s transient dynamics, as given by the solution of 
Eq. § with the nonequilibrium initial condition 

P(0) = k3)(93| (6) 

i.e. with the particle’s probability density initially peaked 
on the right of the potential barrier, in the interval {qb, qc) 
(see Fig. [^. This may be experimentally attained by 
preparing the particle in the ground state of an appro¬ 
priate harmonic well centered at the desired position, and 
then releasing the harmonic potential [ 18 ] , 

Before giving the definition of escape time in the 
present context, we define the population of the lower 
(right side) well as the cumulative population of the three 
DVR states from |q4) to Ige)) that is 
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M=4 

which is a discretized version of the probability of pen¬ 
etration of the wave packet through the barrier m- 
During the transient dynamics the populations of the 
metastable states (|gi) and 192)) reach a maximum. Af¬ 
terwards, by tunneling through the potential barrier, the 
population of the metastable well decays, finally settling 
down to a stationary value dependent on the tempera¬ 
ture. 

We consider large asymmetry of the potential, low tem¬ 
peratures with respect to the barrier height, and damp¬ 
ing regimes ranging from moderate to strong (7 > wq). 
Given the above conditions, the relaxation occurs in the 
incoherent regime, with no oscillations in the popula¬ 
tions. As a consequence we may consider the parti¬ 
cle irreversibly escaped from the metastable state once 
Fright (t) has reached a certain threshold value that we 
set to Fright( t) = 0 . 95 . 

Results. — We observe a nonmonotonic behavior of r 
with respect to both 7 and T . In Fig. [^it is shown the 
presence of a peak in r vs 7, whose height and position 
depend on the temperature. 

A comparison between r and Treiax versus 7 indicates 
that the two quantities exhibit roughly the same behavior 
until the peak in r is reached (see Fig. [2>)- At higher 
7, while Treiax continues to increase monotonically, t has 
a sudden fall off at a critical value 7c, dependent on the 
temperature (for example 7c — 0.98 at T = 0 . 352 ). 

This critical value corresponds to a dynamical regime 
in which the population transfer from the initial state to 
the states of the metastable well is inhibited and there 
is a direct transfer to the states of the lower right well. 
In this regime the probability to find the particle in the 
metastable state is practically zero throughout the en¬ 
tire dynamics. Indeed, while Treiax is the time needed 
for the system to reach the equilibrium in the double 
well potential, the escape time is a relevant quantity for 



FIG. 2: (Color online) Escape time t, in units of for 
the initial condition p(0) = |q 3 )((? 3 | (see Eq. fef). (a) T and 
Treiax US a functiou of both damping 7 and bath’s temperature 
T. (b) T as a function of 7 for different temperatures. Inset - 
escape time vs temperature at fixed values of 7 . The variables 
7 and T are given in units of oio and huo/kB, respectively. 


the transient dynamics, involving the crossing of the po¬ 
tential barrier and the depletion of the metastable well. 
Therefore, our analysis applies to the general problem of 
escape from a metastable well, starting from a nonequi¬ 
librium condition. 

The nonmonotonic behavior of t vs 7 can be inter¬ 
preted as the quantum counterpart of the NFS phe¬ 
nomenon observed in classical systems, and may be called 
quantum noise enhanced stability (QNES). 

Another interesting feature is the presence of a slow 
monotonic increase of t for 7 > 7c, which leads to the 
quantum Zeno effect m- 

The behavior of t vs the temperature is characterized 
by a minimum as ksT approaches the tunneling splitting 
^F'4_3 = Ei — (see Fig. [^. This is the signature of 
the thermally activated tunneling^ an experimentally well 
established phenomenon | 20 j . This is better shown in the 
inset of Fig. [2>. 

Finally we wish to point out that our results are ro¬ 
bust against the variation of the potential asymmetry, 
threshold value and initial conditions (initial DVR states 
within the interval {qb,qc), see Figs. and |^. 
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To exemplify this robustness, we give here the results 
for another potential profile and a different initial condi¬ 
tion (see Fig. [^. 
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FIG. 3: (Color online) (a) Potential profile V(q) (see Eq. 0 ) 
for AU = 2.5hujQ and e = 0.35-\/ Mhcug. Horizontal lines: the 
first 8 energy levels. Vertical lines: position eigenvalues in the 
DVR. The dashed curve is the initial probability distribution 
|v['(a:, 0 )p. Here Aif 7,6 = O.Uhujo, Ai? 6,5 = 0.58?icao, Ai? 5,4 = 
O.lfitJo. (b) Escape time r, in units of for the initial 
condition shown in panel (a). 

The definition of t is the same as for the previous case 
i.e. Pright(T) = 0.95. In this situation however, due to 
the different number of energy levels considered, the right 
well population is defined as Pright{t) = 

The escape time displays qualitatively the same fea¬ 
tures as for the first configuration, even if now the initial 
wave packet is centered close to the top of the potential 
barrier. In particular r has a nonmonotonic behavior as 
a function of both 7 and T. The minimum of r vs T is at 
T ~ 0.27 jkB, which is the average value of the three tun¬ 
neling splittings AE'rfi, and AE^^i (see Fig. [^). 

Moreover, for 7 > 7 c, we observe a monotonic increasing 
behavior of r leading to the quantum Zeno effect. 

Summary. — We have found that in nonequilibrium 
dynamics the escape time from a quantum metastable 
state exhibits a nonmonotonic behavior as a function of 
both the damping, with a maximum (QNES), and the 
temperature, with a minimum at the resonance with the 
tunneling splitting (see Figs. 0 and0. 

We observe stabilization of the quantum metastable 
state due to the dissipation and its interplay with the 


temperature, in the moderate to strong damping regime. 
Moreover, a suppression of the activated escape is ob¬ 
tained by increasing the temperature. The stabilization 
phenomenon associated to our model is within the reach 
of existing experimental technologies such as supercon¬ 
ducting qubits [IS] and optical trapping pij . 

The present model could be used to control the sta¬ 
bility of a trapped particle in atomic optics, precision 
spectroscopy and optical communication. 
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